Asymptotic Expansions of Integrals 

Lectures Fourteen and Fifteen 

In the last lecture, we discuss the method of stationary phase which is applicable to the integral 
(8.36). 

If the integral (8.36) has neither points of stationary phase nor contributions from the endpoints, the 
method of stationary phase fails. As an example, consider calculating the asymptotic form of 

(•GO —lAX 

/(A) = ^ dx. A » 1. 

The function /(A) is the Fourier transform of (1 + x^)"''^. Identifying u{x) with -x, we see that, since 
u'{x) = -1 =^ 0, this integral has no points of stationary phase. Nor does this integral have finite 
endpoints. Therefore, the method of stationary phase is not useful for the problem. In order to find 
the asymptotic form of this integral, we must use the saddle point method whcih will be presetly 
discussed. 

D The Saddle Point Method 

The Laplace method enables us to deal with the integral of the form of (8.8), where the coefficient of 
X in the exponent is real. The method of stationary phase enables us to deal with the integral of the 
form of (8.36), where the coefficient of A in the exponent is purely imaginary. 
In this section, we shall deal with integrals of the form 

I{X) = f e'^^'^h{z)dz, (8.52) 

J c 

where c is a contour in the complex z-plane and where /(z) and h{z) are analytic functions of z. This 
problem is interesting in its own right but, as we have mentioned, it may also help us to evaluate the 
asymptotic form of integrals of the form (8.36) in which there are neither finite endpoints nor points 
of staionary phase. 
Let 

f{z) = u(x,y) + iv{x,y), 
then the exponential function in (8.52) is equal to 

piXJlz) — piXu{x,y) p-Xv{x,y) 
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We see that the coefficient of A in the exponent above has both a real part and an imaginary part. 
In section 8B, we have treated the case in which u is zero. And in Section 8C, we have treated the 
case in which v is zero. We now treat the case with neither u nor v equal to zero. 
The magnitude of the integrand is sharply peaked around the point zo at which v is minimum. But if 
u' is not equal to zero, the integrand also oscillates rapidly in the neighborhood of zo, and the 
resulting cancellation wipes out entirely the contribution of zo- 

Let us deform, by Cauchy's integral theorem, the contour c into another contour c'. If m happens to 
be a constant on the contour c',then 
/(A) = e'^ f e-^''hiz)dz, 

c 

where a is the value of u on c'. Thus there is no cancellation due to the rapid changes of the phase 
factor of the integrand and the dominant contribution to the integral comes from the point zo on c' at 
which V is minimum. The asymptotic form of the integral can be obtained with the use of Laplace's 
method. 

The point zo may be an endpoint of c' . In such a case, we may obtain the leading term of the integral 
by expanding the integrand around the endpoint and carry out the integration over an exponential 
function. 

The point zo may also be an interior point of c' . We shall argue that, if zo is an interior point of c', all 

partial derivatives of u and v must vanish at zo- 

We begin by observing that, since v is minimum at zo along the curve c', the variation of v in the 
tangential direction of c' vanishes at zo- Also, since u is constant on c', the variation of u in the 
tangential direction of c' is zero. The latter means that, by the Cauchy-Riemann equation, the 
variation of v in the direction orthogonal to the tangential direction of of c' is zero. Thus both partial 
derivatives of v vanish at zo- By the Cauchy-Riemann equations, so do the two partial derivatives of 
M at Zo- Indeed, 

/(zo) = 0. 

The point zo is known to be a saddle point for the integral (8.52), as we shall explain further below. 
Since, along the contour c, v is minimum at zo and u is stationary, the dominant contribution to the 
integral over c comes from the saddle point zq. The asjmptotic form of the integral can be obtained 
by expanding the integrand around zo and carrying out the Gaussian integrals. 
As we deform the contour, we may encounter singularities of the integrand. The Cauchy integral 
theorem is then no more apphcable and we may not deform the contour beyond the singularity. Thus 
the singularity may contribute a dominant contribution to the integral. 

In summary, for the integral of the general form (8.52), the possible points of contribution are the 
singular points, the endpoints, and the saddle points. 

There may be more than one saddle points or singular points in the problem. Not all of them 
contribute to the asymptotic form of the integral. We must decide what contour to deform to and 
which saddle points or singular points are relevant. We shall show how to do this below. 
Doing this will require a little understanding of complex analysis. Here we'll add to our previous 
discussions on this topic. 

A curve on which m is a constant is called a level curve of u. A curve on which v is a constant is 
called a level curve of v. We shall show that, with certain notable exceptions, the level curves of u 
intersect orthogonally with the level curves of v. 

Problem for the Reader: 
Consider in the complex plane the function 

Az)=z^ = {x + iy)\ 
Draw the level curves of u as well as those of v for this function. 

Answer 
We have 

u{x,y) = x'^ -y^, and v{x,y) = 2xy. 
We draw the level curves of u and the level curves of v =constant in the figure below. 
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Figure 8.8. 

Question for the Reader: 
Do the level curves of u and those of v intersect orthogonally? 
Answer 

Yes, as we can see from Fig. 8.8, the level curves of u generally intersect orthogonally with those of 

V. 

But there is one notable exception. There are two level curves of u and two level curves of v 
intersecting at the origin. The level curves of u and those of v do not intersect orthogonally at the 
origin. Instead, the two level curves of u intersect orthogonally at the origin. So do the two level 
curves of v. We note also that the derivative of f{z) vanishes at the origin, which is where the 
aforementioned unsual behavior occurs. 

As a general rule, if u and v are harmonic conjugates of each other the level curves of u and those of 
V intersect orthogonally, . To see this, we express the infinitesimal variation of u{x,y) as 

du = Uxdx + Uydy. 

Now on a level curve of u, the value of u does not vary. Therefore we have 

Uxdx + Uydy = 0, 

where dx and d^ are the variations of x and y between two neighboring points on the level curve of u. 
This says that Vw is perpendicular to the tangent vector of the level curve of u. 
Similarly, Vv is perpendicular to the tangent vector of the level curve of v. 
Now 

— > — > 
(Vm) • (Vv) = UxVx + UyVy = 0, 

where the last step is obtained with the use of the Cauchy-Riemann equations (2.9). Therefore, at the 
point a level curve of u and a level curve of v intersect, the tangent vectors of these two curves are 
othogonal to each other. 

But there are exceptions. The points of exception are the points at which the derivative of f{z) ^ 
vanishes. At each of such points, all of the partial derivatives of u and v vanish. Thus both Vm and Vv 
are null vectors at such a point and neither possess a direction. While the scalar product of Vm and W 
remain to vanish at a zero of/, it is due to the vanishing of the gradient vectors, and not due to the 
orthogonality of the gradient vectors. 

Note also that two level curves of u generally do not intersect. Let a level curve be given by 

u(x,y) = c. 

And let {xo,yo) be a point on this curve. Now by the impUcit function theorem, the equation above 
has a unique solution in a neighborhood of {xo,yo) with the form 

y = y(x), 

provided that Uy{xo,yo) * 0. Similarly, it has a unique solution in a neighborhood of (xo,3'o) with the 
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form 

X = xiy), 

provided that Ux{xQ,yo) * 0. Thus, unless both Ux and Uy vanish, there can only be one level curve 
passing through the point (xo, jo)- But at the point zo at which/ (zo) = 0, both Ux and Uy vanish, and 
there is no reason why two level curves of u cannot intersect at zo- Indeed, we see in the example 
above that two level curves of u do intersect at the origin. Also, two level curves of v in this example 
also intersect at the origin. 

We have learned in calculus that if the derivative of a function of a real variable x vanishes at xq, 
then it has an extremum at xq provided that its second derivative does not vanish. This is quite 
different for u{x,y) and v{x,y), which are functions of two variables. Near the point zo at which the 
derivative ofjXz) vanishes, we have 

Az)-Az,) -/'(zo)(z-zo)V2. 

Therefore, if/'(zo) + 0, the behavior of/(z) near zo is the same as the function z^ near the origin. 

Problem for the Reader: 
Let v{x,y) = 2xy, the imaginary part of z^. What is the behavior of v on the line x = y and the line 
X = -y"? What does v look Uke near x = y = 01 

Answer 

The line x = y is a level curve of u. On this line v=2x^, and is minimum at the origin. 
The line x = -y is also a level curve of u. On this line v = -2x^, and is maximum at the origin. 
If we consider v{x,y) as a function of x and y in a three-dimensional space, it resembles the surface 
of a saddle near the point x = y = 0, where the line x = j is in the front-back direction of the saddle, 
and the line x ^ -yis in the sidewise direction of the saddle. The function v{x,y) is neither 
absolutely maximum nor absolutely minimum at the origin. Thus the origin is a saddle point of v. 
Similarly, the origin is a saddle point of m = x'^ -y^. 

If u and V are the real part and the imaginary part of an analytic function the derivative of which 
vanishes at zo, then zo is a saddle point of u and that of v. 

By definition, the value of a function v{x,y) does not vary on any of its level curves. Thus it varies 
most rapidly in the direction orthogonal to its level curve. Since the level curves of u intersect those 
of V orthogonally, v varies most rapidly on a level curve of u, and u varies most rapidly on a level 
curve of v. 

As a consequence, the function e~^^ varies most rapidly on a level curve of u. Now we have seen that 
there are two level curves of u passing through the point zo- We have also seen that the function 
v{x,y) is minimum on one of these curves, and maximum on the other curve. Thus the function 
g-Av(z) decreases most rapidly as one leaves zo on the former curve. This curve is called the path of 
steepest descent. Similarly, e~^^^^^ increases most rapidly as one leaves zo on the latter curve, which 
is called the path of steepest ascent. 

Consider now the integral in the general form of (8.52). If c is a level curve of u(x,y), the phase of 
the exponential function of (8.52) is a constant and can be taken outside of the integral. Also, the 
magnitude ^"''^' (^•y) varies most rapidly on c. Thus the integral is of the form of (8.8), and the 
dominant contributions come from the point at which v is minimum. It is therefore a good idea to 
deform the contour of the integral of (8.52) into a level curve of u, for then the asymptotic form of 
the integral can be obtained with the Laplace method. 

Take the example of (8.47), in which /(z) = z^. We cannot easily go from -1 to 1 on a curve of 
constant u. Let us therefore take the path ACDEFB in the figure. Along the path AC, the phase 
u(x,y) is a constant and v{x,y) is minimum at the endpoint A. Thus the dominant contribution to the 
integral along AC is given by a small neighborhood of width IM near the endpoint A, and can be 
calculated by the Laplace method. The result is the same as we have obtained by using the formula 
(8.39) and keeping only the contribution from the endpoint A. The contour CD is not a level curve of 
u. Instead, it is a level curve of v, on which v = 1. The integrand on CD is as small as e~^, hence the 
integral over CD will be neglected. Next we consider the contour DE. This contour is a path of 
steepest descent with the saddle point at the origin. Thus we may again use the Laplace method to 
evaluate the contribution of the origin, and the answer is the same as that calculated with the use of 
(8.45). The contribution of the contour EF, which is not a level curve of u, will be neglected as the 
integrand is exponentially small on this contour. The dominant contribution of the contour FB, 
which is a level curve of u, comes from the point B. Applying the Laplace method gives us the same 
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result obtained with the use of (8.39). Thus, by deforming the contour of integration, we succeed in 
verifying (8.48) with the apphcation of the Laplace method. 




Figure 8.8. 

We note that it will not be fruitful to deform the contour from the real axis to the curve AC'D'E'F'B. 
For if we do this, we will find that on the curve AC, v is minimum not at A but at C. Thus the 
dominant contributions to the integral over AC come from the point C, where the integrand is of 
the order of e^. This large term is cancelled by an equally large term with the opposite sign coming 
from the integral over CD'. Also, while the contour D'E' passes through the saddle point z = 0, this 
contour is one of the steepest ascent and the contributions of this contour come not from the saddle 
point z = but from the endpoints D' and E'. Indeed, the integral over each of the segments of the 
curve AC'D'E'F'B has a very large leading term. These terms cancel one another and it is extremely 
difficult to determine what remains after such extensive cancellations. 

Therefore, as we start from an endpoint of the contour and try to deform the contour to a level curve 
of M, be sure to go in the direction v increases, not the other way around. 

It is possible to deform the contour into a sum of curves all of which are level curves of u. Instead of 
taking the contour ACDEFB to go from A to B, we may go from A to C on the level curve of u as 
before, but instead of making a turn at C, we follow this level curve all the way to infinity. We then 
change over at infinity to the line x = y, and go from one infinite end of this straightline to another 
infinite end of this straighthne. Then we change over at infinity to the level curve on which F and B 
he to arrive at B. In this way m is a constant on each of the segments of the contours. 
Applying the Laplace method to each of the contours, we may obtain not only the leading term but 
also the next leading terms. To do this, we merely have to keep more terms in the expansions of the 
integrand iiround the points of contribution and carry out the integration. 
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